Deuteron stripping and pickup involving halo nuclei ^^Be and ^^C 
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Three-body calculations of (d,p) and (p, d) reactions involving one-neutron halo nuclei ^^Be and 
^'C are performed using the framework of Faddeev-type scattering equations. Important effects 
of the optical potential nonlocality are found improving the description of the experimental data. 
The obtained values for the neutron spectroscopic factor are consistent with estimations from other 
approaches. 
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I. INTRODUCTION 



Deuteron stripping and pickup reactions {d,p) and 
(p, d) constitute an important tool for extracting nuclear 
structure information such as spectroscopic factors or 
spin/parity. Of special interest are cases involving exotic 
weakly bound systems such as one-neutron halo nuclei 
assumed to consist of a core with a mass number A and 
a neutron n; well known examples are ^^Be, ^^C, and 
-'^^C. The reaction d + A — > p-l- {An) and its time reverse 
p + {An) d-\- A may therefore be described in a three- 
body model {p,n^A). These reactions have been ana- 
lyzed by several approximate methods such as distorted- 
wave Born approximation (DWBA) and various DWBA- 
type adiabatic approaches [11, H, |j, 0, Q or continuum- 
discretized coupled-channels (CDCC) method 
Recently, also the application of exact three-body Fad- 
deev/Alt, Grassberger, and Sandhas (AGS) scattering 
theory [10, to three-body nuclear reactions has be- 
come possible [H, [3] due to a novel implementation 14 1 
of the screening and renormalization method 
for including the long-range Coulomb force between 
charged particles. Faddeev/AGS framework, though be- 
ing technically and computationally most complicated 
and expensive, has an advantage that, once numeri- 
cally well-converged results are obtained, all discrepan- 
cies with the experimental data can be attributed solely 
to the shortcomings of the used optical potentials (OP) or 
to the inadequacy of the three-body model. Recent com- 
parison [l^ of Faddeev/AGS and CDCC results revealed 
that CDCC is indeed a reliable method to calculate d-\-A 
elastic and breakup cross sections but may lack accuracy 
for transfer reactions such as p -f ^^Be ^ d-\- ^°Be. No 
benchmark with adiabatic approaches is performed yet. 
Furthermore, the nonlocal optical potential (NLOP), so 
far, has been included only in the Faddeev/AGS frame- 
work fis'] where an important effect of the OP nonlocality 
was found in deuteron stripping on ^^C and ^^O, improv- 
ing the description of the experimental data. Therefore 
we expect the OP nonlocality to be significant also in 
transfer reactions on one-neutron halo nuclei which we 
aim to study in the present work using momentum-space 
AGS equations. 

The theoretical framework is shortly recalled in Sec. [Til 



The results for (d,p) and (p, d) reactions involving ^^Be 
and ""^^C nuclei for which the experimental data is avail- 
able are presented in Sec. IIIII The summary is given in 
SecHVl 



II. AGS EQUATIONS 



The AGS equations [ll| 

Upo.{Z) = ~5po. G^\Z) + ~5p^T^{Z) Go{Z)Ujc.{Z), 

7 

(1) 

are a system of Faddeev-like coupled integral equations 
for the transition operators Uf3a{Z); their on-shell ma- 
trix elements {il'f3\Uf3a{E + iO)\il!a) at the available three- 
particle energy E are amplitudes for all scattering pro- 
cesses (elastic, inelastic, transfer, breakup) allowed by 
the chosen Hamiltonian H = Hn + y^., v. 
the free Hamiltonian and v-y 
7 in odd-man-out notation. In Eq. ([1]) Spa = 1 — <^/3aj 
Go{Z) = {Z — Hq)^^ is the free resolvent, and 



Ho + J2j ^1 where Hq is 
the potential for the pair 



T^{Z) ~ + v^Go{Z)T^{Z) 



(2) 



is the two-particle transition matrix. The channel states 
I ■07) for 7 = 1, 2, 3 are the eigenstates of the correspond- 



ing channel Hamiltonian Hj = Hq 



with the energy 



eigenvalue E; thus, lipj) is a product of the bound state 
wave function for pair 7 and a plane wave with fixed on- 
shell momentum q.^ corresponding to the relative motion 
of particle 7 and pair 7 in the initial or final state. 

The AGS equations ([T]) are applicable only to 
short-range interactions. Nevertheless, the long-range 
Coulomb force between charged particles can be included 
in this framework usin g t he method of screening and 
renormalization p^. [TgI. TlTj where one has to solve AGS 
equations with nuclear plus screened Coulomb potential 
to obtain the Coulomb-distorted short-range part of the 
transition amplitude; the convergence of the results with 
the screening radius has to be established. The method 
has been successfully appHed to proton-deuteron elastic 
scattering and breakup [13, [3| and to direct nuclear 
reactions dominated by three-body degrees of freedom 

MM. 
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We use momentum-space partial-wave basis. The 
screened Coulomb potential and most of the nuclear po- 
tentials are given in coordinate space; they have to be 
transformed to the momentum space where local and 
nonlocal potentials are treated in the same manner. The 
AGS equations then become a system of integral equa- 
tions with two continuous variables which are the abso- 
lute values of Jacobi momenta. The technique of numer- 
ical solution is described in Refs. [iqI. f20l. [2l| . 



III. RESULTS 

The dynamic input to AGS equations are the poten- 
tials v-y for the three pairs of particles. For the np in- 
teraction we take the realistic CD Bonn potential [2^ . 
In order to describe reactions involving one-neutron halo 
nucleus we need a real nA potential that reproduces the 
bound state spectrum of nucleus (An). This is a stan- 
dard choice in p+ (An) reactions where the pA potential 
is complex and is taken at the proton lab energy. As dis- 
cussed in Ref. il8(], d+A ^p+{An) and p+{An) d+A 
reactions are related by time reversal provided the energy 
in the center of mass (cm.) system is the same. There- 
fore we can calculate the latter one using the standard 
Hamiltonian with the nucleus {An) being in its ground 
or excited state and apply the detailed balance to obtain 
the observables for the former one. This is equivalent to 
calculating the d + A —^p + (An) transfer with a real nA 
potential and complex pA interaction which is a nonstan- 
dard choice. Nevertheless, it provides quite a reasonable 
description of the d+A elastic scattering, as demon- 
strated in Ref. [l^ and in the present work. Therefore 
we use a real partial-wave dependent nA potential that 
has local central and spin-orbit parts, 

Vy{r) = -VJ{r, R,a)+(T-L V,o -^f{r, R, a), (3) 

r dr 

with /(r,i?,a) = [l + exp((r-i?)/a)]~i and i? = roA^/^. 
The nA^Be potential is taken from Ref. [23| while the one 
for n-^'^C is constructed in the present work assuming 
standard values for vq ~ 1.25 fm, a — 0.65 fm, and Vso — 
5.5 MeV and adjusting Vc to the binding energy of ^^C 
ground and excited state and to the neutron separation 
energy of ^*C. For completeness the parameters of both 
potentials are given in Table U while the resulting values 
for binding energies are listed in Table [Til Pauli forbidden 
bound states are projected out as described in Ref. p^ . 

For the hadronic part of the pA interaction we take 
the NLOP of Giannini et al. [l^l which is based on Wat- 
son's multiple scattering theory; the nonlocality arises 
mainly due to the fully off-shell two-nucleon transition 
matrix [2^, [2^ . In the configuration space the NLOP is 
parametrized as 



^;^(r',r) = H,{x)U,{y) + Ht{x)Ut{y) 
+ Hs{x)Vs{y)(T -l. 



A-2Z 
A 



(4) 



tons 



with X = [r' - r|, y = \v' + r|/2, Hj{x) = 
exp {—x^/Pj), and Z being the number of pro- 
in the nucleus A. The values of the nonlocality 
parameters are /3c = 1.015 fm, (3t — 1.633 fm, jSs = 0.789 
fm. The nonlocality of the OP mainly Both Uc{y) and 
Ut{y) contain real volume and imaginary surface parts, 
while the spin-orbit part Vs{y) is real; all of them are 
parametrized in the standard way using Woods-Saxon 
functions and their derivatives [2J|, i.e., 

Uj{y) = - V,f{y,aR) - AtW,f{y,ai)[l - f{y,ai)], 



Vsiy) 

f{y,ak) 



Vs 



2 df{y,aR) 



V dy 

[i + c^v{{y-RN)/ak)]-^ 



(5a) 
(5b) 
(5c) 



with Rn = lAQA^/^ fm, an = 0.57 fm, a/ = 0.54 + 
0.0032^1 fm, Vc = 85 MeV, Vt = 127 +UZA-^/^ MeV, 
Vs = 9.1 MeV, Wt = 13 MeV, and Wc = wn[1 - 
exp(-0.05£;)] MeV where E = + 1.08- 1.35^^-^/3 
and is the proton energy (MeV) in the cm. frame. 

We adjust wn to improve the description of the experi- 
mental pA scattering data in the energy regime relevant 
for the considered three-body reactions. In the case of 
p-^^C we fit the NLOP to the data at proton lab en- 
ergy Ep = 14.5 MeV We need the p-^°Be poten- 
tial in a broader energy r ang e than Ep = 12 — 16 MeV 
where experimental data [28| is available. Since the lo- 
cal energy-dependent OP by Watson et al. [l^ provides 
rather satisfactory description of that data as demon- 
strated in Ref. [28| and confirmed by our own calcula- 
tions, we adjust the NLOP to the predictions of Watson 
OP. The obtained values for wn are give n in Table IIIIl 
other parameters are taken from Ref. [24| . This time we 
do not use the local OP obtained from NLOP by equiv- 
alence transformation [23|; instead, the Watson OP that 



TABLE I: Partial-wave dependent strengths of central and 
spin-orbit parts of n-^'^Be and n-^*C potentials, all in units 
of MeV. Other parameters: ro = 1.25 fm, while a — 0.67 and 
0.65 fm for n-^"Be and n-^*C , respectively. 

L K(n-^°Be) V;o(n-^°Be) Vc(n-^^C) Ko(n-'^C) 





odd 

even, > 2 



56.413 
42.498 
56.413 



11.953 
5.38 



50.29 
46.13 
49.14 



5.50 
5.50 



TABLE 11: Binding energies (MeV) of the bound states cor- 
responding to the potential parameters of Table HI Pauli for- 
bidden bound states that are removed are marked with *. 



Is 



1/2 



2si/2 



lP3/2 lPl/2 



Ids 



5/2 



"Be 28.730* 
^^C 28.194* 



0.503 6.812* 0.183 

1.218 11.522* 8.175* 0.479 
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TABLE III: The parameter uijv of NLOP adjusted to the p- 
data and p-^"Be Watson OP predictions at given proton 
lab energies Ep (all in units of MeV) . 



Ep Wiv 

p-^°Be 9.0 40.0 

p-^°Be 12.0 36.0 

p-^°Be 20.9 27.0 

p-"Be 35.3 18.0 

p-^*C 14.0 28.0 



is often used also in DWBA-type and CDCC calculations 
serves as a local reference potential. 

Although some parameters of those pA potentials vary 
with energy, in the calculations of present paper based 
on the Hamiltonian theory for three-body scattering they 
are taken at a fixed energy, corresponding to the proton 
lab energy in the p + (An) reaction. Therefore the cho- 
sen Hamiltonian with complex pA potential prevents the 
calculation oi d + A ^ n+ (Ap) and p + {An) ^ n + {Ap) 
reactions since the nucleus (Ap) is not bound. These 
reactions can be described [30'| allowing for an explic- 
itly energy-dependent pA potential which becomes real 
at negative relative pA energies and reproduces bound 
states of the nucleus (Ap). However, such an approach 
is not free of theoretical complications due to the ab- 
sence of the Hamiltonian theory as discussed in Ref . [3^ ; 
we therefore refrain from using it in the present work. 
Furthermore, the results of Ref. [1^ indicate that the 
effect of the energy-dependence in the pA potential on 
the (d, p) and {p, d) cross section is rather insignificant 
at small cm. scattering angles up to Qc.m. ~ 40°, al- 
though it may become sizable for large angles where the 
cross section itself is small. 

The interaction between np, nA, and pA pairs is in- 
cluded in partial waves with pair orbital angular momen- 
tum L < 3, 6, and 12, respectively, and the total angular 
momentum is J < 25; the results are well converged. 
The pA channel is more demanding than the nA channel 
due to the screened Coulomb force, where the screening 
radius i? ~ 10 to 12 fm for the short-range part of the 
scattering amplitude is sufficient for convergence. 

Examples for the cross section in two-body p + ^"Be 
and p + ^^C and three-body d + ^'^Be and d + ^^C elastic 
scattering are presented in Figs. [T] and [2] and compared 
with the available experimental data. Two-body results 
indicate that adjusting only one parameter wj^ of NLOP 
is sufficient to achieve the quality in fitting the pA data 
comparable with the one of the Watson OP. However, the 
three-body d + A data seems to be reproduced better by 
the NLOP, although it also fails for d + ^°Be at large an- 
gles. There is no p -I- (An) elastic scattering data in the 
considered energy region; however, as demonstrated in 
Ref. [l^, the theoretical predictions ior p+A and p+ {An) 
elastic scattering observables are strongly correlated and 
therefore also the nonlocality effect on p + (An) elastic 
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FIG. 1: (Color online) Differential cross section divided by 
Rutherford cross section for p-l-^^Be elastic scattering at Ep = 
12 MeV and d + ^°Be elastic scattering at E^. = 12 MeV. 
Predictions of NLOP (solid curve) and Watson OP (dashed- 
dotted curve) are compared with the experimental data from 
Ref. [SI. 
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FIG. 2: (Color online) Differential cross section divided by 
Rutherford cross section for p-|- ^*C elastic scattering at Ep = 
14.5 MeV and d + "C elastic scattering at Ed = 14 MeV. 
Curves as in Fig. [1] and data from Refs. [27l, [3H. 



cross section is extremely small. Thus, the observed dis- 
crepancy in p-f ^^Be elastic cross section at higher energy 
[Hi [13] cannot be resolved by the NLOP. 

Results for d + ^"Be — > p + ^^Be transfer at deuteron 
lab energy Eij = 12 and 25 MeV are presented in Fig. [3] 
for the ^^Be ground state 1/2+ and in Fig. 2] for the 
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FIG. 3: (Color online) Differential cross section for d+^"Be —>■ 
p+ ^^Be transfer to the '^^Be ground state 1/2"'' at Ed — 12 
and 25 MeV as function of the cm. scattering angle. Curves 
as in Fig. [1] and the data are from Refs. [2^. [33|. 



spectroscopic factors below 1. This is believed to be due 
to ne glec ting the ^°Be core excitation in the ^^Bc bound 
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FIG. 4: (Color online) Same as Fig. O but for transfer to the 
^^Be excited state l/2~. 



excited state 1/2 . The corresponding proton lab en- 
ergy in the time-reverse reaction p + ^^Be ^ d + ^'^Be 
is Ep = 9.0 (8.7) and 20.9 (20.5) MeV for the ground 
(excited) state, respectively. At least for the transfer to 
the ground state 1/2+ the NLOP reproduces the shape 
of the data clearly better than Watson OP. In all cases 
theory overpredicts the data and thereby yields neutron 



state [3J|. For the ground state 1/2+ , the spectroscopic 
factor value around 0.45 at 12 MeV is in agreement with 
the DWBA-type result of Ref. [|, while at 25 MeV the 
value around 0.75 is consistent with most other estima- 
tions collected in Ref. [s^. Slightly larger values of the 
neutron spectroscopic factor are obtained for the excited 
state 1/2", around 0.6 and 0.85 at 12 and 25 MeV, re- 
spectively. In Fig.[n]results for p-|-^^Be d-l-^^Be trans- 
fer at Ep ~ 35.3 MeV are presented. Again significant 
difference between NLOP and Watson OP predictions 
is found. The values for the ^^Be ground state neutron 
spectroscopic factor obtained using the NLOP from (d, p) 
reaction at Ed = 25 MeV and {p, d) reaction at Ep = 35.3 
MeV are consistent with each other, in contrast to those 
obtained using the Watson OP. 

Results for d+^'^C ^p + ^^C transfer at Ed = 14, 16, 
and 17 MeV are presented in Fig. [5] for the ^^C ground 
state 1/2+ and in Fig. [7] for the excited state 5/2+. The 
corresponding proton lab energy in the time-reverse re- 
action ^ d-h"C is Ep = 12.0 (11.2), 13.9 (13.1), 
and 14.8 (14.0) MeV for the ground (excited) state, re- 
spectively. The theoretical predictions vary smoothly 
with the energy while the data from Ref. [s^ at 16 
MeV are sig nificantly lower than other two sets from 
Refs.j3ll[36{ at 14 and 17 MeV as already pointed out in 
Ref. [5|. This raises serious doubts on the normalization 
of the 16-MeV data. In all cases NLOP reproduces the 
shape of the data better than Watson OP, especially for 
the transfer to the ground state 1/2+ in the region of the 
first minimum and second maximum. At 14 MeV and 
at 17 MeV up to cm. scattering angle 6c. m. « 60° also 
the quantitative description of the data by the NLOP is 
quite satisfactory. The extracted neutron spectroscopic 
factor for the ^^C ground state 1/2+ is close to 1, in 
agreement with the DWBA-type result of Ref. ^ , while 
for the excited state 5/2+ it is slightly below 1. 
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FIG. 6: (Color online) Differential cross section for d + ^''C p + ^"'C transfer to the ^^C ground state 1/2+ at Ed = 14, 16, 
and 17 MeV. Curves as in Fig. [1] and the experimental data are from Refs. [3ll. Issl. [3^ . 
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FIG. 7: (Color online) Same as Fig. [G] but for transfer to the ^^C excited state 5/2^ 



In addition we performed the nearside-farside decom- 
position of the cross section [s^ • Although at small an- 
gles where the cross section is largest the nearside usu- 
ally dominates, the nonlocality effect may be equally im- 
portant for both nearside and farside cross sections and 
therefore the observed nonlocality effect is a result of a 
complicated interference between nearside and farside. 

The results of the present work together with those of 
Ref. [13 indicate that the OP nonlocality effect on trans- 
fer observables depends quite strongly on the reaction en- 
ergy and internal orbital angular momentum of the (An) 
nucleus but far less on its binding energy and mass. The 
dependence on the internal orbital angular momentum 
is evident, and the dependence on the reaction energy is 
seen most clearly by comparing d+ ^"Bc p+ ^^Be and 
p + "Be ^ d + ^°Be results at = 12 and 25 MeV 
and Ep — 35.3 MeV given in Figs. [3] and [5l However, 
as shown in Figs. [3] and [51 the OP nonlocality effect is 
very similar in {d,p) reactions at Ed = 12 and 14 MeV 



leading to the ground state 1/2+ of ^^Be and ^^C, re- 
spectively. Furthermore, the OP nonlocality effect in the 
p + "Be ^ d+ i^Be transfer at Ep = 35.3 MeV is very 
similar to the one shown in Figs. 5 and 7 of Ref. [l^ for 
{d^p) reactions at comparable cm. energies leading to 
the excited state 1/2+ of ^'^C and ^^O nuclei with the 
binding energies of 1.857 and 3.272 MeV, respectively. 
There is also a mutual similarity between OP nonlocal- 
ity effects given in Figs. 5 and 6 of Ref. [ll| for {d,p) 
reactions around Ed = 30 MeV leading to 5/2+ state of 
^^C and ^^O nuclei; the corresponding binding energies 
are 1.092 and 4.143 MeV. 

As in Ref. [l^ we checked that the OP nonlocality ef- 
fect observed in transfer reactions is insensitive to small 
variations of the parameter wm and therefore is not a 
consequence of only approximate on-shell equivalence be- 
tween NLOP and Watson OP. We also studied sensitivity 
of the results with respect to np and nA interactions. We 
performed additional test calculations using local AV18 
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potential [33| for np] compared to mildly nonlocal CD 
Bonn potential the differences in the predictions were 
found to be entirely negligible. Varying parameters tq 
and a of the nA potential within reasonable limits, ±0.05 
fm, or, in the case of n-^^Be, using i-independent Vso, 
and adjusting the strengths Vc and Vso to reproduce the 
same binding energies, leads to the differences in trans- 
fer results that are still considerably smaller than the 
observed OP nonlocality effect. Thus, the cross section 
in the considered transfer reactions depends mostly on 
the pA interaction. 

IV. SUMMARY 

We performed calculations of (d, p) and (p, d) transfer 
reactions involving one-neutron halo nuclei ^^Be and ^'^C. 
Exact three-body scattering equations in the AGS form 
were solved and the Coulomb interaction was included 
using the method of screening and renormalization; well 
converged results were obtained. With respect to dynam- 
ics, they are mostly sensitive to the proton-core interac- 
tion. For the first time in reactions with halo nuclei a 



nonlocal optical potential was used and important non- 
locality effects were found. The OP nonlocality effect on 
transfer observables depends mostly on the reaction en- 
ergy and internal orbital angular momentum of the halo 
nucleus and is rather insensitive to its binding energy and 
mass. The NLOP accounts for the experimental trans- 
fer data better than the local Watson OP and provides 
a more consistent description over a broader energy and 
angular range, especially for the 1/2+ initial/final state 
of the halo nucleus. The values for the neutron spectro- 
scopic factor obtained using the NLOP agree quite well 
with estimates based on other approaches, in contrast to 
those obtained using the nearly on-shell equivalent Wat- 
son OP. 
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